Basically, how bad of an approximation can we do to the interpolation and still get away with 10% accuracy to the grid points? 10K, 0.05 dex in logg and [Fe/H] ?
As a comparison we need a high quality spectrum generated using our normal interpolation methods.
The SPEX tests are different than the TRES tests. This is because we are going from R=500,000 to R=2000, and thus the convolve first then interpolate will bring out more errors. The interesting thing, though, is that the error from doing convolve-interpolate is a different spectrum than the interpolation envelope. This is good. It means that this might be simply like adding interpolation noise (that the GP should pick up) rather than actually adding in any specific bias.
Towards the 4000 Ang example, we might have to be careful about interpolating after convolution. In this case, it might be better to pre-sample the grid to higher resolution w/ the splines, then convolve.
The areas where the HA is much lower than the LA are generally also continuum regions. This is good news, too.
In [1]:
%matplotlib
import matplotlib.pyplot as plt
import numpy as np
from StellarSpectra.spectrum import DataSpectrum
from StellarSpectra.grid_tools import SPEX, HDF5Interface
from StellarSpectra.model import Model
from grid_accuracy import AccuracyComparison
from IPython.display import display
In [2]:
myDataSpectrum = DataSpectrum.open("../../data/Gl51/Gl51.hdf5")
myInstrument = SPEX()
In [3]:
myComp = AccuracyComparison(myDataSpectrum, myInstrument, "../../libraries/PHOENIX_submaster_M.hdf5",
"../../libraries/PHOENIX_SPEX_2300.hdf5",
{"temp":3000, "logg":4.5, "Z":-0.5}, {"temp":30, "logg":0.1, "Z": 0.1})
fig = myComp.plot_quality()
display(fig)
In [11]:
myComp = AccuracyComparison(myDataSpectrum, myInstrument, "../../libraries/PHOENIX_submaster_M.hdf5",
"../../libraries/PHOENIX_SPEX_2300.hdf5",
{"temp":3000, "logg":4.5, "Z":-0.5}, {"temp":10, "logg":0.05, "Z": 0.05})
fig = myComp.plot_quality()
display(fig)
In [12]:
myComp = AccuracyComparison(myDataSpectrum, myInstrument, "../../libraries/PHOENIX_submaster_M.hdf5",
"../../libraries/PHOENIX_SPEX_2300.hdf5",
{"temp":3000, "logg":4.5, "Z":-0.5}, {"temp":1, "logg":0.01, "Z": 0.01})
fig = myComp.plot_quality()
display(fig)
In [13]:
myComp = AccuracyComparison(myDataSpectrum, myInstrument, "../../libraries/PHOENIX_submaster_M.hdf5",
"../../libraries/PHOENIX_SPEX_2300.hdf5",
{"temp":4000, "logg":4.5, "Z":-0.5}, {"temp":30, "logg":0.1, "Z": 0.1})
fig = myComp.plot_quality()
display(fig)
In [14]:
myComp = AccuracyComparison(myDataSpectrum, myInstrument, "../../libraries/PHOENIX_submaster_M.hdf5",
"../../libraries/PHOENIX_SPEX_2300.hdf5",
{"temp":4000, "logg":4.5, "Z":-0.5}, {"temp":10, "logg":0.05, "Z": 0.05})
fig = myComp.plot_quality()
display(fig)
In [15]:
myComp = AccuracyComparison(myDataSpectrum, myInstrument, "../../libraries/PHOENIX_submaster_M.hdf5",
"../../libraries/PHOENIX_SPEX_2300.hdf5",
{"temp":4000, "logg":4.5, "Z":-0.5}, {"temp":1, "logg":0.01, "Z": 0.01})
fig = myComp.plot_quality()
display(fig)
In [ ]: